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Abstract 



The asymptotic behaviour of a family of inhomogeneous scalar 
field cosmologies with exponential potential V = Ae k< ^ is studied. By 
introducing new variables we can perform an almost complete analysis 
of the evolution of these cosmologies. Unlike the homogeneous case 
(Bianchi type solutions), when k 2 < 2 the models do not isotropize 
due to the presence of the inhomogeneities. 
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1 Introduction 



In the last years there has been a renewed interest in studying scalar field 
cosmological solutions mainly due to the prominent role that scalar fields 
would play in the early stages of the universe driving a period of inflationary 
expansion. Scalar fields can be thought as arising from the the low-energy 
limit of string theory or as coming from the dimensional reduction in Kaluza- 
Klein models. 

One of the most interesting questions that has been analysed is the so- 
called cosmic no-hair conjecture which states that the dynamics of a scalar 
field with a realistic potential could isotropize the space time smoothing 
out the inital anisotropies and inhomogeneities, giving at late times a FRW 
universe. Initially, the cosmic no-hair conjecture was proved for space times 
with a cosmological constant [[!]]. Since during the "slow- rolling" epoch the 
scalar field behaves as a cosmological constant, that result was supposed to 
be applicable to scalar field cosmologies. However, it was soon realised that 
the picture would change when the dynamics of the scalar field is taken into 
account 0. The study of the evolution of cosmological models with scalar 
field has been done mainly for homogenous solutions (Bianchi type) and, in 
particular, when the potential of the field is of Liouville type (exponential 
potential) ||. Exponential potentials arise in dimensional reduction theories 
and it is well known that they lead to power law inflation in FRW models. 

Although there have been several attempts to study the evolution of the 
scalar field in inhomogeneous models @], there is a lack of results concern- 
ing how spatial inhomogeneities affect the dynamics of the scalar field. It 
was shown for a particular family of exact inhomogeneous solutions with ex- 
ponential potential that there are cases for which the scalar field does not 
guarantee that the models inflate or isotropize ||. Since a generic relativistic 
solution near the initial singularity is neither homogeneous nor isotropic it is 
worth analysing the effect of spatial inhomogeneities on the isotropization. 

In a previous paper (hereafter Paper I) M we started the study of the 
asymptotic behaviour of a particular class of inhomogeneous metrics de- 
scribed by a diagonal Gi line-element. In that work we analysed a particular 
family of G<i self-similar solutions with exponential potential. The main con- 
clusions were that the asymptotic behaviour of the model does not depend 
on the parameters describing the spatial inhomogenity but rather on the pa- 
rameter k, associated with the slope of the exponential potential and the 
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behaviour, essentially, was the same as that found for homogeneous solu- 
tions, i.e: when k 2 < 2 the models isotropized. It was surmised that these 
conclusions could be explained in terms of the "weak" spatial inhomogeneity 
and that the introduction of more "strong" inhomogeneities could result in a 
different behaviour with respect to the homogeneous case. Since our analysis 
is based only on the study of the equilibrium points, the conclusions we draw 
have to be considered at the conjecture level, until a more complete analysis 
of the structure of the field equations for G2 metrics is done. 

In this paper we extend the analysis of Paper I to a more general class 
of G2 exponential potential cosmologies and we show how in those cases 
for which the homogeneous solutions isotropized (when k 2 < 2), here the 
inhomogeneities prevent the isotropization. Since our analysis is based only 
on the study of the equilibrium points, the conclusions we draw have to 
be considered at the conjecture level, until a more complete analysis of the 
structure of the field equations for G 2 metrics is achieved. 

It is known that scalar field homogeneous cosmological solutions with an 
exponential potential evolve towards self-similar solutions and it was conjec- 
tured that the same is true for G 2 solutions J7|. In this paper we prove that 
the asymptotic metrics of the class of solutions studied here are self-similar. 

The plan of the paper is the following: in Section 2 we will describe the 
metric and the assumptions we make in order to simplify the equations. A 
new set of variables is introduced in such a way that the phase-space of 
the system is bounded for most of the cases considered. In Section 3 the 
equilibrium points of the system and the correponding solutions are studied. 
Finally the conclusions are explained in Section 4. 

2 The metric and the dynamical systems 

In this paper we shall be considering spatially inhomogeneous G2 cosmologies 
admitting a two-parameter abelian group of isometries for which the line 
element is written as: 

ds 2 = e F (-dt 2 + dz 2 ) + G (e p dx 2 + e~ p dy 2 ) , (1) 

with all functions appearing in the metric depending on both variables t and 
z. In paper I we considered that the transitivity surface function G was a 
homogeneous function, depending only on t, and that the rest of the metric 
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functions were separable. In the present case we would like to extend that 
analysis allowing the function G to depend on the variable z as well, keeping 
the separability of the metric functions. 

The local behaviour of the models described by metrics of the form ([!]) 
is determined by the gradient of the so-called transitivity surface function 
G |. A globally null or spacelike gradient of this function corresponds to 
plane or cylindrical gravitational waves. If the gradient is globally timelike, 
however, the metric describes a cosmological model with spacelike singu- 
larities and if the character of the gradient changes throughout space time 
it describes colliding gravitational waves or cosmologies with timelike and 
spacelike singularities. 

As the matter source we shall consider a minimally coupled scalar field 
with exponential potential: 

V{(j>) = Ke k(t> . A > (2) 

For vanishing k this will be equivalent to a cosmological constant term. The 
corresponding stress-energy tensor will read: 



'ab = 9,a<P,b - 9ab 



(3) 



It is well known || that this stress-energy tensor is equivalent to a perfect 
fluid as long as the gradient of the scalar field is timelike. 

For the model under consideration the resulting Einstein and Klein-Gordon 
field equations are: 
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-^' 2 = -0 2 + 0' 2 . (9) 

We want to find the general behaviour of the models under consideration 
and, in particular, their behaviour near the initial singularity and at late 
times. The treatment of the asymptotic behaviour of the inhomogeneous 
metrics was initiated in [|?| by extending the method developed for homoge- 
neous Bianchi type solutions with a perfect fluid [ lOj] . This method is based 



on the introduction of a new set of adimensional variables constructed from 
the fluid quantities. For inhomogeneous metrics the new variables satisfy a 
system of partial differential equations and the equilibrium points are con- 
sidered as dynamical equilibrium points. Since the source of our metrics, 
the scalar field, is not globally equivalent to a perfect fluid we can not apply 
this method and we have to tackle the problem of resolving the asymptotic 
behaviour of system (4)-(9) in a different way. In paper I we analised a par- 
ticular family of solutions of the above system. In particular, we considered 
a family of solutions for which the inhomogeneity was linear and introduced 
new dynamical variables in such a way that the phase-space was compact 
and, therefore, we could describe completely the asymptotic behaviour of 
the system. We will follow the same approach as that already used in pa- 
per I. By fixing the spatial dependence of the metric functions we will write 
this system of partial differential equations as a set of ODE's and we will 
use standard dynamical system techniques to find out the structure of the 
phase-space. 

To obtain the desired set of ODEs we have made a number of simplifying 
assumptions on the metric functions: 

i) As in |TTJ we write the scalar field in terms of a new function ip as: 



= ~log(G)+V>- (10) 

Substituting (10) into the Klein-Gordon equation (4) we find that ip satisfies 
the same equation as function p. This allows us to choose the function ip 
proportional to p: 

ip(t, z) = mp(t, z) (11) 
where m is an arbitrary constant. 
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ii) to simplify the equations further we consider separability of the metric 
functions and in particular we assume the ansatz, already used by different 
authors @,[|nj: 



G(t,z)=T{t)Z(z) (12) 
F(t,z) = f(t) + f 1 (z) (13) 
eP(M) = Q(t)Z(z) n . (14) 

These are ad hoc assumptions, made with the only motivation of simpli- 
fying the Einstein equations. They have been suggested by previous works 
on inhomogeneous metrics with perfect fluid source. 

With these assumptions, the two functions depending on the spatial vari- 
able z, that is, fi(z) and Z(z) are inmediately fixed by the equations. In 
fact, we see that Eq(6) separates into: 

Z" TQQ (Q\ 2 2 

where nea 2 stands for an arbitrary separation constant. Resolving the above 
equation we get: 

{Acosh(az) + Bsmh(az) e = 1 
Az + B e = (16) 

A cos(az) + B sin(az) e = — 1 

where A and B are arbitrary constants of integration. 

From Eq.(|5|) when A ^ or from Eq(^) when A = we find that: 

fi{z) = 2a\ogZ(z) + logL>. (17) 

When A = 0, a and D are arbitrary constants, but when A ^ 0, D is still an 
arbitrary constant but 

2a = C = k 2 /2 - kmn. (18) 

The rest of the equations separate in a straightforward way giving rise to 
new separable constants and a few relations betwen them. We shall not give 
here the details of the calculations but only summarize them. 
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When Z(z) 7^ e ±az we obtain from Eq.(8) or from Eq.(9) the following 
constraint on the parameters: 

2« = n 2 M 2 + C-^ + ^, (19) 

where we have defined M 2 = 1/2 + m 2 . When A = the above equation 
fixes the value of a, but when A 7^ it becomes, from (18), 

n 2 M 2 = I + ^! (20) 

2 4 K ' 

The last relation allows us to write m in terms of n and k. The two solutions 
are: 

yp + 2 -2n 2 
m ± = ± (21) 



which also gives C: 



k 



C± = -(k^Vk 2 + 2-2n 2 ). (22) 



From now on we assume that Eqs(19)-(22) hold. That means that, in 
the exponential case, we are only considering those solutions for which the 
constants verify the above relations. In the general exponential case all the 
solutions are self-similar and the constant m is not restricted by relation (|2l|) . 
Since this case will not add any new insights to the main conclusions of the 
paper we leave the corresponding analysis to be discussed elsewhere. We are 
also assuming that A^O. Solutions with A = shall be limiting cases from 
which all the solutions with A ^ start their evolution. 

Finally the line element and scalar field have the form: 

ds 2 = De f Z 2a (-dt 2 + dz 2 ) + TQZ^ 1+n) dx 2 + ^-Z^dy 2 (23) 
—k 

4> = — logT + mlogQ + [ran - k/2) \ogZ. (24) 

Where the function Z(z) is defined in (16). 

Among the solutions described by the line element (23) we can identify 
the following known solutions, all of them correponding to the case when 
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a = 0: If Z{z) = sin(az), or cos(a;z), De* = QT and n 2 = 1, then the 
solutions reduce to the Kantowski-Sachs models; if Z(z) = e ±az , then one 
obtains that: when n 2 = (23) reduces to a type Bianchi V solution, when 
n 2 = 1 to a type Bianchi III and when n 2 ^ 0, 1 to type Bianchi VI. Finally, 
if Z(z) =constant (23) is a type Bianchi I solution. 

We can see also from Eq.(24) that when mn = k/2 the scalar field is 
homogeneous and, therefore, can be cast in a perfect fluid form. In this case 
the metric (23) belongs to the classes of solutions studied by different authors 



Another important class of solutions belonging to the line element (23) 
is that obtained when k = m = (from Eq.(20) n 2 = 1). In this case the 
scalar field vanishes and the source of the metric is a cosmo logical constant. 

Once the spatial dependence has been fixed, we have to find the time de- 
pendence of the metric that is determined by the three independent functions 
T(t), Q{t) and fit). We introduce, for later convenience, a new function h(t) 
defined as: 

h(t) = f + kmlogQ. (25) 
The Einstein equations, with the above assumptions, reduce to: 

^ = ea 2 + 2 AD exp (h - y log t) , (26) 

Q 2 TQ (Q\ 2 

Q =nta -TQ + \q) ' (2?) 

(28) 

(29) 



k 2 

= 2ADexp(h - — logT) + 2ea 2 . (30) 

From the original system of PDE's we have obtained a system of five 
ODE's. These equations separate into three dynamical equations and two 
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constraint equations. Since equation (|26|), however, depends on T and its 
second derivative the associated dynamical system will have four first order 
equations and two constraint equations. On the other hand, the constraint 
fl3"0|) will serve to eliminate the parameter ea 2 from the system. 

To put this system of equations in the usual dynamical system form, 
we make an appropriate change of variables. We want the new variables 
to be normalized in such a way that they remain bounded near the initial 
singularity. 

From Eq.(||) we see that an appropriate set of new variables (3, $ and T 

is: 



with 



M- 



b-f + dh 
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b-f + dh 



T = V2AD 



exp (± 
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f logT 



bj; + dh 



b 2 
b = 



_ |2~fc 2 | 
4 ' 

d=l 



d 2 



\2-k 2 \ 



k 2 ^2 
k 2 = 2 



(31) 



(32) 



T is related with the potential of the scalar field, it vanishes when the 
scalar field is massless, (3 is related with the time dependent part of the 
function p in the metric ([!]) and $ depends on /. Unlike the Bianchi type 
solutions, where the kinematical quantities of the source are used as vari- 
ables, here, since the scalar field is not globally timelike, we can not use 
that approach. The reason for defining the variables as in (|31"D is to get a 
bounded phase space (at least for a large subset of the solutions). Thus there 
is no particular physical meaning for these variables, although another set of 
variables with some physical meaning could perhaps be found. 

In terms of the new variables (31) and substituting ea 2 with the aid of 
Eq.(30), Eq.(26) becomes: 



T 
T 



^b' 



1 b $ 



Y 2 -f3 2 



8 V 



(33) 



To decouple this equation from the rest of equations of the system we define 
a new time r, through the following expression: 



T 

b— + dh = — 
T 1 



T 



$ T 



(34) 
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From the general behaviour near the initial singularity found by Belinskii 
et al [13], we can assume that our models, for constant z, behave near the 
initial singularity like a Kasner model, in the sense that when t — > 



T ~ t, f ~ logQ ~ ~ logi. (35) 

This means that f ~ 1/t when t — > which guarantees the regular behaviour 
of the new variables near the initial singularity. To see how r behaves far 
from the singularity we write Eq. flSOf ) in the following way: 



k 2 

f 2 = 2ADexp(h- — logT) + 2ea 2 + M 2 

+ d2h2 +( b2 + k ^-)(f) 2 (36) 

when k 2 ^ 2 and 

f 2 = 4A J Dexp(/i-logT)+4ea 2 + 2M 2 (^j + h 2 + (^j (37) 
when k 2 = 2. 

As long as e = 0, 1 all the terms in the right hand side of the above two 
expressions are positive and, therefore, f is different from zero for all values 
of t. Because of this we can restrict our attention to f > which implies 
that our new time r will grow monotonically with t. Since f ~ 1/t for t ~ 0, 
the range of r goes from — oo to +oo. On the other hand, the situation with 
e = — 1 is completely different. Since there is nothing preventing f from 
being zero at some finite time t, there will be solutions whose trajectories in 
the phase-space of the new variables will escape to infinity in a finite time. 
That means that, in this case, r will no longer serve to monitor the evolution 
in the coordinate time t. 

We are now in a position to write the equations in terms of these variables. 
Since the change of variables does not depend on k in a continuous way, 
the system of equations does not have the same form for all the values of 
k. First we write the constraint equations (29) and (30) that in these new 
variables take particularly simple expressions. The first constraint equation 
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(29) relates $ and (3: 



A-k 2 



k 2 = 2 



$ = 



2 + yW^4f3 
k 2 



k 2 >2. 



As for the second constraint Eq.(30) it takes the simple form: 

2ea 2 



L(p, $, r) = i - p 2 - $ 2 - r 2 + H(p, $, r) 



where 



# = k 2 <2 

H = <f>-1 k 2 = 2 

\2 1.2 



(38) 



(39) 



(40) 



^ H = -2($ - l) 2 A; 2 > 2. 

Finally, substituting (31) in Eqs.(26)-(28) we get the dynamical equations 
for the different ranges of the parameter k 



k 2 <2: 



r' 



i 

46 

r 

46 



(3 - 3$ - 26 2 ) L + 2 [l - (26 2 + l) $ 
-3L + 2 (26 2 + l) (l - T 2 ) - 2$ 



(41) 



k 2 = 2: 



/3' 

r' 



(1 - 4/3) L - 2/3r 2 



= - (3 - 4$) L + (1 - 2$) P 



!' |2L-- + 2r 2 



(42) 
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k 2 >2: 



r' 



^- ([(3 + Ab 2 ) (1 - $) - 26 2 ] L + 2 [l - (l + 26 2 ) $] T 2 ) 



(43) 



46 



(3 + 46 2 ) L - 2 (l - 26 2 ) + 2 (l - 46 2 ) $ + 2 (l + 26 2 ) 



The prime stands for the derivative with respect to the new time r and L is 
defined in (39). 

From (39)-(40) we see that there are three different invariant subspaces, 
one for each value of the parameter e. The first is the closed surface L = 
which corresponds to the linear case e = 0. The other two subspaces are 
separated by the surface L — 0. The invariant subspace for e = 1 is a 
compact region limited by the surfaces L = and T = and, finally, the 
invariant subspace for e = — 1 is the unbounded region outside L — 0. 

The most satisfactory scenario to give the explicit form of the asymptotic 
behaviour of the models is the one in which we have variables which are 
bounded for all times. This is guaranteed for both e = and e = 1 since we 
have seen that in these two cases f remains always positive and the phase- 
space is bounded. As explained above this is not the case for e = — 1. We 
will see from the numerical calculation of the trajectories that, as expected, 
for e = —1 some trajectories go to infinity in a finite time. Probably one 
can find a different set of variables for which this part of the phase space 
becomes bounded, although it is unlikely that this can be done for the whole 
phase-space. 



3 The equilibrium points and the associated 
metrics 

We are interested in looking for results concerning no-hair theorems, thus 
for each of these dynamical systems we will look for the equilibrium points 
and their stability. This is important since attractors in phase-space will be 
related to the behaviour of the models in the limit r — > 00 (t — > 00) while 
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sources, the unstable fixed points, will be related to the behaviour near the 
initial singularity r — > — oo (t — > 0). 

Before giving the explicit expressions of the equilibriun points of the 
system, it is interesting to investigate whether the equilibrium points rep- 
resent self-similar solutions. It has been shown that the homogeneous but 



anisotropic models (Bianchi models) evolve towards self-similar solutions [10 
and it has been conjectured that the same behaviour is shared by the G 2 met- 
rics . We will prove that the solutions considered in this paper do approach 
self-similar solutions. Here is a sketch of the proof of this statement. 

When e = 0, 1 the equilibrium points are characterized by finite constant 
values of the variables j3, $ and T, then for an equilibriun point equation 

(33) takes the form: 

f = (1 - W) (pj (44) 

where W is a function both of the constant values of j3, $ and T and the 
parameters of the model. This equation has only two types of solutions 
depending on whether or not W is zero. When W 7^ the solution is of the 
form T = t l / w while when W = it is an exponential, that is, T = e wt , 
where W is another constant. For the power law solution T = t x l w , equations 

(34) and (39) force e = which in turn forces Z(z) to be linear in z. 
Writing down the equations for a homothetic vector £, 

£^9ab = 2g a b, (45) 

we find that for both solutions, exponential and power law, there is at least 
one homothetic vector with components £ a = (do, a±x, a,2y, 0) in the exponen- 
tial case and £ a = (b t, bix, b 2 y, b z) in the power law case, where both a m and 
b m are constants (the details of this calculation are given in the Appendix) 
To give a qualitative behaviour of the solutions of the system (41)- (43) 
it is important to obtain the equilibrium points and their local stability. Let 
us first note that the system admits the symmetry T — > —T. Therefore, the 
study of the equilibrium points will be restricted to V > 0. They can be 
found explicitely and are given by: 
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Pi and P 2 



p 2 + $ 2 = 1 r = k 2 < 2 

^2 + $ 2_ $ = o r = k 2 = 2 

/3 2 + (1 - $) (1 - 3$) = T = k 2 >2 

These points lie on the surface L = corresponding, therefore, to e = 0. 
All of them are the sources of the entire phase-space, i.e: the attractors 
at early times. This set of points reduces to two separate points when the 
constraint equation (38) is taken into account. These are the two points given 
by the intersection in T = of the surface L = and the plane defined by 
the first constraint equation (38). For k 2 > 6 one of them becomes a saddle, 
attracting the trajectories lying on L = 0. Since T = 0, the corresponding 
metrics describe massless scalar field cosmologies. The line element and the 
scalar field are: 

ds 2 = Dt c ^E^(-dt 2 + dz 2 )+t 1+c *E] +n dx 2 + t 1 - c *E l - n dy 2 

where 

E z = Az + B (47) 
C^l + f (1 + C ± ) 
C 2 = n(l ±v /=5). 

The — and + signs appearing in the definition of C<i correspond to P 1 and P2 
respectively. When n 2 = 1 and mn — k/2 C± = and it is easy to see by a 
simple coordinate transformation that (46) is a LRS Bianchi type I solution. 

Pa 



= Wp! $ = 1 _ | 6 2 r = k 2 < 2 
P=\ $ = § r = 0A; 2 = 2 

o _ 2bvrw (t) _ 3+2fe 2 r — n P -> 9 
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This point corresponds, again, to a massless scalar field but wiht e = 1. 
It is a saddle for all k, attracting the trajectories with e = 1 in T = 0. The 
metric and the scalar field corresponding to this point are: 



ds 2 = De ^+C±)at E C ± (_ dt 2 + dz 2 ) + ( e at E J 1+n d% 2 + ^at^ — dy 2 



1-n 



= 



(at + log E z 



(48) 



where 



£7* = A cosh a,2 + B sinh 



and A and B are arbitrary constants. 

By making the change: T = e at E z , Z = e at E' z when A > B we see that 
(47) is an anisotropic Bianchi type I space time. (For A < B the model 
has spacelike curvature singularities). In the special case when n 2 = 1 and 
mn — k/2 this turns out to be Minkowski flat metric. 



Pa 



/3 = $ = 1 



■n _ 2fev / I+& 2 " 
1+262 



k 2 <2 



1+26 2 



r _ 26y / I^F 7,2 v. o 
1 _ 1+26 2 K > Z 



For this point e = 0. It does not exist for k 2 = 2 or k 2 > 6. For k 2 < 2 
this is the attractor for trajectories with e = and e = 1 (also for some of 
the trajectories with e = —1). For 2 < k 2 < 6 it is a saddle, attracting only 
the trajectory with e = 0. In this case: 



ds 2 = 

= 



k 2 



2k 



-.t 



4 



2A 



-dt 2 + dz 2 ) + El +n dx 2 + E L ~ n dy 



2-k 2 



log 



k 2 -2 
2 V / 6^P' 



^logE 2 , £ z = As + £. (49) 



When n 2 = 1 and mn = k/2 (C± = 0) it is easy to see that this is a flat 
FRW space time. 
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Pk 



P = 7^m $ = l-26 2 T = - 7 2L= k 2 <2 



1 r = ± 

2 2 
o _ 2b 3 VTW (fs _ l+2fc 2 r _ 2fcy / T+b 2 " L2 ^ o 



/3 = $ = \ r = ± k 2 = 2 



l+56 2 +46 4 ^ 1+4&2 - 1 - l+56 2 +46 4 

The solution corresponding to this point has e = — 1 when A; 2 < 2, e = 
when k 2 = 2 and e = 1 when k 2 > 2. It is a saddle for k 2 < 2 and a sink for 
A; 2 > 2. When k 2 > 2 this is the attractor for the trajectories with e = 1 and 
when fc 2 = 2 is the attractor for both trajectories with e = and e = 1. The 
metric and the scalar field are: 

ds 2 = A2 - €a2 e C 1 t E C ±{ _ dt 2 + dz 2 ) 

2A 2 v ' 

+ e [(A 2 +nea 2 )/A]t E l+n dx 2 + g [(A 2 -nea?) / 'A\t gi-n dy 2 

1 (C 1 t + C±logE z ) (50) 



k 



where 



k 2 A 2 T ea 2 k^k 2 + 2 - 2n 2 2 A; 2 -2„ 2 
Cl = 2A ' = e ~k 2 +~2 

Acosh(az) + i?sinh(a2;) e = 1 
E, = { Az + B e = 

Acos(az) + Ssin(az) e = — 1 

where A, A and are arbitrary constants. 

This solution was previously found in [O . As in the previous cases, when 



n 2 = 1 and mn = k/2 it is homogeneous: when k 2 < 2, it is a Kantowski- 
Sachs space time; when fc 2 = 2 it is a flat FRW model while for k 2 > 2 it is 
a LRS Bianchi type III solution. This last case can be seen by making the 
coordinate transformation: 

X = -dx-lnE z , Z = -^^e~ dx , 

d l E z 

where d is a constant related with the constants appearing in the solution. 

To complete the description of the evolution of the solutions we have 
solved the system (41)-(43) with the constraint equation (38) by numerical 
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integration. The phase-space diagrams so obtained are plotted in Fig.l for 
different ranges of the parameter k. The bold line represents the invariant 
subspace L = 0. This line splits the phase-space into two separate invariant 
subspaces: the inner region correponds to e = 1 (L > 0) and the outer 
region to e = — f (L < 0). The two fixed points Pi and P2 are located in 
which constitutes the common points of the two invariant subspaces defined 
by L = and T = (subspace of massless scalar field solutions). P 3 , lies in 
the invariant subspace r = 0, L > 0. P 4 lies in the line L = and exists 
only when k 2 < 2 and when 2 < k 2 < 6. It aproaches the point Pi as k 2 — > 6 
and the point P5 as k 2 — * 2. Finally, P5 does not stay in the same invariant 
subspace for all k. As a matter of fact, this point is in the region defined by 
L < when k 2 < 2, is on L = when k 2 = 2, and lies in L > when k 2 > 2. 

4 Conclusions 

In this work we have considered a class of inhomogeneous solutions with a 
minimally coupled scalar field with an exponential potential. This class is 
more general than that studied in paper I. Therefore, we can analyse the 
effect produced by complex initial inhomogeneities on the dynamics of scalar 
field cosmologies and specifically on the isotropization of the models. 

When the inhomogeneity was linear, as that described in paper I, or when 
the solutions were homogeneous (Bianchi type), the dynamics of the solutions 
was completely determined by the value of the parameter k and it was shown 
that for k 2 < 2 all the solutions reached homogeneity and isotropy but when 
k 2 > 2 while all the models still attained homogeneity only a subclass of 
measure zero reached isotropy. 

We have found the equilibrium points of the dynamical systems associated 
with part of the solutions (23) and we have shown that all of them are self- 
similar. Except the point P 3 , they correspond to inhomogeneous space times 
with a scalar field configuration that cannot be related to a perfect fluid 
source. However, in the special case when the scalar field is homogeneous 
(n 2 = 1 and mn = k/2), they are all homogeneous perfect fluid space times. 
The fixed point P3 represents a homogeneous perfect fluid space time for 
every value of the parameters defining the particular model. 

We can conclude that although the dynamics is still determined by k, as in 
the former cases, the result refering to the isotropization or homogeneization 
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of the space time is completely different. Actually, when k 2 < 2 solutions 
with e = 0, 1 evolve towards the point P4 which is inhomogeneous. Although 
the point P 3 is an attractor for a subset of solutions, it is a saddle point, 
attracting only the massless scalar field solutions. Hence, the inhomogeneities 
prevent in this case the solutions from isotropizing. Only in the particular 
case when n 2 = 1 and mn = k/2 (corresponding to a homogeneous scalar 
field) we do have a similar situation to the previous cases. When k 2 > 2 the 
point P5 becomes the attractor for the solutions. 

In the region where e = — 1 the situation is quite different since the vari- 
ables used are not bounded. If k 2 < 2 the stable manifold of the fixed point 
P5 splits this region into two separate regions. The solutions with trajecto- 
ries taking place in the inner region are attracted by the inhomogeneous fixed 
point P 4 , whereas those on the stable manifold approach the inhomogeneous 
solution P 5 . Those with trajectories in the unbounded outer region escape 
to infinity. When k 2 > 2 this separation does no longer exist and all the 
solutions with (e = —1) escape to infinity. 

Points Pi and P2 describe the solutions "close" to the initial singularity 
t = 0. On both points the potential vanishes since near the singularity the 
kinematical term is the leading one in the stress-energy tensor. By making 
a coordinate transformation the line element and the scalar field can be cast 
into the form: 

ds 2 = -a (z)dt 2 + ai (z)t 2pi dx 2 + a 2 (z)t 2p2 dy 2 + a 3 (z)t 2p3 dz 2 , 
(f> = q\nt + b(z), 

where pi and q are constants verifying 

Pi + P2 + P3 = vX + ?i + ?i + Q 2 = 1- 



This agrees with [14] in the sense that the above metric could represent the 
asymptotic form of a "general" scalar field cosmological solution near the 
singularity. 

It is important to note that solutions corresponding to a cosmological 
constant can be analysed by letting k = m = 0. In this case all the solutions 
evolve to homogeneous space times. 

As we mentioned in section 2, in order not to extend this paper the general 
exponential case will be considered in a separate paper. It is relevant just to 
note that the asymptotic behaviour in that case is similar to that found here 
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in the sense that even when k 2 < 2 the equilibrium points of the system are 
inhomogeneous. 
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A Appendix 

In this appendix we are going to show that all the fixed points of the system 
under study are indeed self-similar solutions. By inspecting the equilibrium 
points listed in Section 3 it is easy to see that for the fixed points the functions 
and Q(t) are powers of the function T(t): 

e m=T(t) Cl Q{t)=T{t) c \ (A.l) 

where C\ and C 2 are constants depending on the equilibrium point we are 
considering. 

— # 

The equations for the existence of a homothetic vector £, taking into 
account Eq.(A.l), reduce to the following equations: 



c^e + 2a^e + 2d e = 2 


(A.2) 


(l + C 2 )^° + (l + n)^ 3 + 2^ 1 = 2 


(A.3) 


(i - c 2 )^e + (i - n)^e + 2d 2 e = 2 


(A.4) 




(A.5) 


DT Cl Z 2a d x e - T 1+C2 Z 1+n d ^ = 


(A.6) 
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DT Cl Z 2a d 2 £° - T l -°' 2 Z l - n d^ 2 = (A.7) 
d 3 e - d o e = (A.8) 
T 2C *Z 2n d 2 £ 1 + d^ 2 = (A.9) 

DT Cl Z 2a d^ 3 + T 1+c ' 2 Z 1+n d 3 ^ = (A. 10) 

DT Cl Z 2a d 2 £ 3 + T 1 ~ C2 Z 1 ~ n d 3 £ 2 = (A. 11) 



We have shown that the function T(t) in the equilibrium points verifies 
the Eq.(44) giving two different solutions: 
r)W = 0=>lf = k = constant. 

Looking at the equations (A.2)-(A.ll) one can try the following solution: 

e = constant ? = £\x) £ 2 = £ 2 (y) f = 0- (A.12) 

Equations (A.5)-(A.ll) are automatically satisfied and equations (A. 2)- 
(A.4) determine completely £: 

- f 2 d - C 2 - 1 d + C 2 - 1 1 

e = I dife' — d — x ' — d — y ' °; • (A - 13) 

II) W + =► % = ^ and f = \ 

In this case equations (A. 2) and (A. 5) suggest a solution of the form: 

e = a ?=z\x) e = e(y) e = cz, (a.u) 

where c is a constant to be determined. 

As before, the actual form is completely determined by equations (A. 2)- 
(A.4), the remaining equations being automatically satisfied. Then we ob- 
tain: 



£ = 



w 



d + W(C 3 + 2) 
d + d - 1 



2t, 



d — d 



W 



+ 1 + d 



n x, 



W 



+ l + C 3 + n)y,2z 



(A.15) 
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Figure Caption 



Figure 1: Phase-space diagrams are drawn for (a) k 2 = 1, (b) k 2 = 2, 
(c) k 2 = 4 and (d) k 2 = 10. Each diagram represents the typical form of the 
phase-space for the four significant ranges of the parameter k. T is shown on 
the vertical axis and f3 on the horizontal axis. The bold line stands for the 
separatrix L = 0, the arrows show the direction of increasing time, and the 
dots represent the fixed points, which are named as in the body of the paper. 
We can easily distinguish the four important invariant subspaces: T = 
(massless scalar field), e = — 1 (the unbounded region outside L = 0), e = 
(the separatrix L = itself) and e = 1 (the bounded region iside L — 0). 




(c) 2< k < 6 



(d) k > 6 



